LYAPUNOV FUNCTIONALS FOR THE ENSKOG EQUATION 



ZHENGLU JIANG 



Abstract. Two Lyapunov functionals are presented for the Enskog equation. One is to de- 
scribe interactions between particles with various velocities and another is to measure the 
distance between two classical solutions. The former yields the time-asymptotic convergence 
of global classical solutions to the collision free motion while the latter is applied into the 
verification of the stability of global classical solutions. 



1. Introduction 

We are concerned with two Lyapunov functionals to describe the time-asymptotic behaviour 
and the stabihty of global classical solutions to the Enskog equation without any external 
force for a hard sphere gas. In the absence of any external force, the Enskog equation is as 
follows: 

where / = f{t,x,v) is a one-particle distribution function that denpends on the time t G R+, 
the position x G and the velocity v £ R^, and Q is the Enskog collision operator whose form 
will be addressed below. Here and throughout this paper, R+ represents the positive side of the 
real axis including its origin and R^ denotes a three-dimensional Euclidean space. 

The collision operator Q is expressed by the difference between the gain and loss terms 
respectively defined by 

Q~^{f){t,x,v) = a'^ / f{t,x,v')f{t,x — auj,w')B{v — w,uj)du!dw, (1-2) 

Q^{f){t,x,v) = / f{t,x,v)f{t,x + auj,w)B{v — w,uj)dujdw (1-3) 

In equations JESI) and Sl = {uj £ : uj{v-w) > 0} is a subset of a unit sphere surface S 
in R^, a is a diameter of hard sphere (a > 0), c<j is a unit vector along the line passing through 
the centers of the spheres at their interaction, {v' , w') are velocities after collision of two particles 
have precollisional velocities {v,w) and B{v — w,uj) = max(0, (v — w)uj) is the collision kernel. 

This equation (|1.H) is indeed the commonly known Enskog-Boltzmann equation and it is 
a modification of the original work of Enskog [5]. There are other different versions of the 
Enskog equation in order that they formally satisfy some such properties as entropy bound and 
consistence with irreversible thermodynamics (see ^ , (HI , ^Hj ) ^Ej ) • It is worth mentioning that 
equation ()1.1() still obeys the conservation laws of mass, momentum and energy |12j . 

As for the Boltzmann equation, two colliding particles obey the conservation laws of both 
kinetic momentum and energy, as follows. 

y^W = y'^U)', V"^ + W'^ = V''^ + W''^ . (1.4) 
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This results in their velocity relations 

v' = V — [{v — ■w)u)]lj, w' = w + [{v — w)uj]uj, (1-5) 

where w G 5^. This implies that B{v' — w', —uj) = B{v — w,io). 

The Boltzmann equation models dilute gases successfully but it is no longer suitable for gases 
with high-density effects. The Enskog equation, a partial differential integral equation of the 
hyperbolic type, is a model first proposed in 1922 by Enskog 5 as the generalization of the 
Boltzmann equation describing the dynamical behavior of the density of a moderately dense 
or high-density gas. There are many results about the global existence and uniqueness of the 
solutions to the initial value problem for the Boltzmann and the Enskog equations without 
external forces. Concerning the Boltzmann equation in the absence of external forces, a global 
solution existence is obtained by DiPerna & Lions |2j for the large data but one cannot yet 
know whether the solution to the problem is unique or not; however, a global existence and 
uniqueness result is shown by Illner & Shinbrot [7j about the solutions to the initial value 
problem to the Boltzmann equation for small initial data in the infinite vacuum. Concerning 
the Enskog equation without external forces, a global existence and uniqueness proof is given 
by Polewczak for near-vaccum data and another one shown by Arkeryd for the large data. 
Some different existence results are also given by Cercignani (see j2], jSl)- With an analysis of 
the well-posedness of the initial value problem in unbounded domains, some global existence and 
uniqueness theorems is obtained by Toscani & Bellomo about the solutions to the Enskog 
equation in the absence of external forces for small initial data with suitable decay to zero at 
infinity in the phase space, and the asymptotic stability of the solutions and the influence of the 
external field have been discussed. For the Enskog equation in the absence of external forces, 
the stability of solutions is first given by Cercignani and the time-asymptotic behaviour 
of solutions in the weighted L°° is then provided by Polewczak (see [Hj , fill j ) • Many other 
results about this subject can be found in the references in the papers mentioned above. 

Two so-called Lyapunov functionals recently constructed by Ha 6^ mathematically neither 
yield the time-asymptotic behaviour in the norm nor recover the stability for global 
classical solutions to the Enskog equation without external forces (see Appendix^. Now there 
are not yet any Lyapunov functionals to restore the the time-asymptotic behaviour and the 
stability of solutions to the Enskog equation in the absence of external forces. The aim of 
this paper is to build two Lyapunov functionals for the Enskog equation. One is to describe 
interactions between particles with various velocities and another is to measure the distance 
between two classical solutions. The former yields the time-asymptotic convergence of global 
classical solutions to the collision free motion while the latter is applied into the verification of 
the stability of global classical solutions. 

The rest of this paper is arranged as follows. In Section ^ some properties of the collision 
operator Q of the Enskog equation are introduced including both the entropy identity and 
the nonincreasing property of the entropy functional. Then in Section |21 a Lyapunov functional 
is constructed and its time-decay property is given together with an application of the time- 
asymptotic behaviour of any nonnegative solution to the Enskog equation (jl.lj) . A different 
Lyapunov functional is defined and its time-decay property, together an application of the 
stability of solution to the Enskog equation (jLlj) . is provided in Sectional A counterexample 
of inequalities given by Ha is finally shown in Appendix 1X1 

2. Preliminaries 

In this section some properties of the collision operator Q of the Enskog equation 1)1. 1|) are 
introduced including both the entropy identity and the nonincreasing property of the entropy 
functional, and the conservation law of mass is shown for the Enskog equation. 
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We consider the Enskog equation with (|1.2|) and Notice that B{v — w,oj) = 

max(0, {v — w)uj). Then we have the fohowing lemma. 

Lemma 2.1. Suppose that Q is a collisional operator as defined by and Let ip{x,v) 

be a measurable function on Y{? >i Y{? . If f ^ Co(R+ x R-^ x R^), then 

/ 'ip{x,v)Q{f)dvdx = — / B{v — w,uj)f{t,x,v)f{t,x + auj,w) 

X [tpi^x, v') — tpix, v) + ip{x + atj, w') — ip{x + auj, w)]dujdwdvdx (2-1) 

and 

/ '4}{x,v)Q{f)dvdx = — I B{v — w,uj)f{t,x — auj,w')f{t,x,v') 

JR3xR3 2 Jr3xR3xR3xS^ 

X [ijj{x — ato, w) — ip{x — ato, w') + tpix, v) — '4>{x, v')]diodwdvdx. (2-2) 



Proof. Denote Ig and // by 
— 

and 



Ig = / B{v — w,uj)ip{x,v)f{t,x,v')f{t,x — auj,w')dujdwdvdx (2-3) 



R3xR3xR3x5' 



■2 



Ii = / B{v — w,u!)il){x,v)f{t,x,v)f{t,x + au),w)dujdwdvdx, (2.4) 

respectively. By (|1.2j) and (|1.3|1 . we have 

il>{x,v)Q{f)dvdx = Ig-Ii. (2.5) 



/R»xR3 

Let us first consider the loss integral Ii . By exchanging v and w and replacing cj with —uj in the 
integral on the right side of (|2.4jl . we can get 

Ii = / B{v — w,u>)il^{x,w)f{t,x,w)f{t,x — aui,v)dujdwdvdx. (2.6) 

JrSxRSxRSxS^ 

Replacing x with x + aw in the integral on the right side of (|2.6() gives 



Ii = a / B{v — w,oj)tp{x + aijo,w)f{t,x,v)f{t,x + au:,w)dujdwdvdx. (2.7) 

JR3xR3xR3x52 

Combining 1)2. 4|) and 1)2. 7|) . we have 
/■ 

Ii = — B{v — w,uj)[ip{x,v)+i/j{x + au;,w)]f{t,x,v)f{t,x + auj,w)du)dwdvdx. 

2 yR3xR3xR3x52 

(2.8) 

Then the gain integral is below considered. By using the properties that B{v' — w' , —lo) = 
B{v — w,ui) and that dwdv = dw'dv' , (|2.3|) can be rewritten as 



Ig = a- 



/ B{v' — w' ,—oj)ip{x,v)f{t,x,v')f{t,x — aLO,w')dujdw'dv'dx. (2.9) 

JR3xR3xR3x5i 



By exchanging v and v', w and w' , and replacing uj with —uj in the integral on the right side of 
(f2T9|) . we have 

Ig = / B{v — w,uj)'4j{x,v')f{t,x,v)f{t.,x + auj,w)dujdwdvdx. (2.10) 

JR3xR3xR3x5i 
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By repeating the analysis of the loss integral Ii discussed above, (|2.1U|) can be rechanged as 
h ~ '^^ / B{v — w^Ljj)il}{x + auj^w')f{t,x^v)f{t^x-\-auj^w)dijjdwdvdx. (2.11) 

JR3xR3xR3xS2 

Combining (|2.1()j) and (|2.1H) . we can get 

Ig = — / B{v — w,uj)['ip{x,v') + ip{x + auj,w')]f{t,x,v)f{t,x + auj,w)dujdwdvdx. 

2 yR3xR3xR3x52 

(2.12) 

Inserting (fTH|) and (|TT^ into (|2ISI) gives (EU). 

We below prove (|2.2|) . By replacing x with x — au in the integral on the right side of (|2.1j) . 
we know that 

ip{x,v)Q{f)dvdx = — / B{v — 'w,u;)f(t,x — auj,v)f{t,x,w) 

R3xR3 2 Jr3xR3xR3xS2 

X [V'(x — aw, v') — ip{x — auj, v) + Tp{x, w') — ij{x, 'w)]dLodwdvdx. (2.13) 

Exchanging v and w and replacing uj with —lo in the integral on the right side of H2.13() . we can 
obtain 

/ 'ip{x,v)Q{f)dvdx = — / B{v — w,uj)f{t,x + auj,w)f{t,x,v) 

JR3xR3 2 Jr3xR3xR3xS2 

X [ip{x + aw, w') — ip{x + au;, w) + ^(x, u') — tpix, v)]dujdwdvdx. (2.14) 
Notice that B{v — w,uj) = B{v' — w', —uj) and that dwdv = dw'dv' . Then (|2.14j) gives 

a^ f 

Tp{x,v)Q{f)dvdx = — / B{v' — w' ,—u))f{t,x + auj,'w)f(t,x,v) 

R3xR3 2 Jr3xR3xR3x52 

X [tp^x + au, w') — ^/'(x + aw, u)) + V'(3;, v') — Tp{x, v)\dujdw' dv' dx. (2.15) 

By exchanging v and v' , w and w' , and replacing w with — w in the integral on the right side of 
(f?T^ . (EH follows. □ 

We can also know that if ijj{x,v) = 1, v, f^/2 in Lemma l2 . 1 1 then 

ip{x,v)Q{f)dvdx = 0. 

R3xR3 

This formally implies the conservation laws of mass, momentum and energy hold for the Enskog 
equation in the absence of external forces, i.e., 

d f ( v'^\ 

1, w, — /(t, X, v)dvdx = 0. 



dt 



R3xR3 



We now introduce a notation /"^(t, x, v) = f{t, x + vt, v). Then the Enskog equation (|l.lj) can 
be also written as 

^%^ = <3(/)*«,:.,.). (2,16) 
Then let us consider two functionals Hsit) and I{t) defined as follows: 

HB{t)= f{t,x,v)lnf{t,x,v)dxdv (2.17) 

JR3xR3 
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and 

f 

I{t) = — B{v—w,LL!)[f{t,x,v)f{t,x+atLi,w) — f{t,x,v)f{t,x — atLi,w)]dujdwdxdv. 

2 Jr3xR3xR3x5;^ 

(2.18) 

They have the following property. 

Lemma 2.2 ( 10 ). Assume that H{t) = Hsit) + fQl(s)ds where Hsit) and I{t) are given by 
\2.11\j and i2.18\) . If f = f{t,x,v) is a nonnegative classical solution to the Enskog equation 
JO) with and then < 0. 

Proof. By H2.17|l . we have 

dHBjt ) _ f , i„ ,#u ^ ..,,9f*{t,x,v) 

dt 

Inserting (|2.1()j) into (|2.19j) gives 
dHB{t) 



[ [l + lnf*{t,x,v)]^^-^^^^^dxdv. (2.19) 



[ [l + \n f*{t,x,v)]Q{f)*{t,x,v)dxdv. (2.20) 

JR3xR3 



dt 

By integrating over the variable x on the right side of (|2.2()j) . it follows that 
dHBjt) 

dt JR3xR3 

By Lemma l2.H (|2.21j) can be rechanged as 



/ {\^\^f(t,x,v)\Q{f){t,x,v)dxdv. (2.21) 

JR3xR3 



dt 



X In 



B{v — w, uj)f{t, X, v)f{t, X + auj, w) 
RSxR^xRaxS^ 

f{t,x,v')f{t,x + auj,w') 



diodwdxdv. (2.22) 

J[t,x,v)f{t,x + au;,w) 

Equation (|2.22j) is also called entropy indentity. Since Inx < x — 1 for x > 0, the estimation of 
the integral on the right side of the entropy indentity l|2.22|) reads 

— — < — / B{v-w,uj) 

dt 2. Jr3xR3xR3x52 

X [/(t, X, v')f{t, X + aw, w') — f{t, X, v)f{t, X + am, w)]du)dwdxdv, (2.23) 
i.e., '^^^l''''^ < ~I{t)- Therefore < ~I{t) +-^(0 — 0- The proof of this lemma is finished. □ 

3. Asymptotic Behaviour 

In this section two new functionals are built. One of them is a Lyapunov functional and its 
time-decay property is related to another functional. This time-decay property can be applied 
into our description of the time-asymptotic behaviour of any nonnegative solution to the Enskog 
equation JHH) with (fO|) and (fO)) . 

Let us begin with two new functionals D = S+ -|- 2)^ and 5 = + d , where 2)^ and 5^ 
are defined as follows: 



^^{f]{t)= [ f*{t,x,v) 

JR3xR3xR3 

/ f^it, X + (v — w)t ^ aoj + T- — ,w)dTduj 

jR+xSi {V-W)UJ 



dwdvdx, (3-1) 
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\{v — x, v)f'^{t, x + {v — w)t =F aw, w)dujdwdvdx. (3.2) 



R3xR3xR3x52 



Obviously, by and we know that J^z^-^s Q^{f)*{t,x,v)dvdx < a^d^[f]{t). For any 

nonnegative solution / to the Enskog equation, D[f]{t) describes interactions between particles 
with various velocities. It can be below found that the time-decay property of D[f]{t) for any 
nonnegative solution / to the Enskog equation leads directly to the time-asymptotic behaviour 
of this solution in L^(R^ x R^) and so T> is called Lyapunov functional. To estimate the time 
decay of D[f]{t), we first have to show the following lemma. 

Lemma 3.1. Let F{x) he an integrahle function on Yi? and v a vector in Yi? . Assume that 
5^ = {wl^o; > 0,u; G S*^} where S"^ is a unit sphere surface in R^. Then 



/ F(x)dx = aM F[auj + r— 
Jr3 JRx5i 



)dTduj 



(3.3) 



where a 7^ 0. 



Proof. Denote x by x = (xi,X2,X3). Using a transformation x = auj + tv/{vuj) where to = 
(sin cos sin sin cos 0), 9 is an angle between v and uj, < 9 < ■|,0<(^< 27r. we can 

know that the Jacobian determinant of this transformation is ~ a^sin^. Lemma l3.ll 

therefore holds. □ 



Combining 1)3. 1|) and Lemma [3.11 we can easily deduce that S)^[/](t) < ^||/||^i for any 
nonnegative integrable function / = f(t,x,v). Furthermore, we also obtain the following time- 
decay property of D [f] (t) . 

Theorem 3.2. Let D and ^ be defined by i'^.l]) and iV. Assume that f = f{t,x,v) is 
a nonnegative classical solution to the Enskog equation / li. ij) with il.S]) and through an 

initial datum /o = fo{x,v) and that f € L^(R? x R'^). then 

d^[f]{t) 



dt 



+ (1-4||/o||li )?[/](*) <o 



(3.4) 



for any t £ R+. 

Proof. We first calculate }P • To do this, we have to deduce that 



/*(t, X, v)f*{t, X + {v - w)t - auj + T- 



V — w 



w) 



dr 



{v — vu)oj 

{y — w)ujf'^(t, X, v)f'^{t, X + {v — vu)t — acj + T 
+Q{f)*{t, X, v)f*{t, x + {v-w)t-auj + T 
+f*{t, X, v)Q{f)*{t, x + {v-w)t-auj + T 



{v — w)uj 
V — w 



,w) 



{v — w)u! 

V — w 



w) 
w) 



(3.5) 



{v — w)oj 

for any fixed variables (x,f, 1/7,0;). In fact, (|3.5j) results from (|2.16|) and the two identities as 
follows: 



dt 



/*(t, X + {v — w)t ~ auj + T 
{v — w)dxf'^{t, X + (f — w)t — auj + T 



V — w 



w 



{v — w)oj 

V — w 



{v — w)uj 



w 
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+Q{f)*{t,x + {v-w)t-auj + T-^ — ^,w) (3.6) 

[V — W)UJ 



and 

(v — w)dxf'^(t, X + (v — w)t — aiv + T- 7 — , 1 

[v — w)uj 

= (v — w)ujdr X + {v — w)t — auj + t- — 

[ (v — W)UJ 

Denote two functionals 3f (i = 1, 2) by 

:ff[/]W= / Q{f)*{t,x,v) 

JRSxRSxRSxS^ xR+ 



(3.7) 



x/*(t, X + {v — w)t =p au; + r- — ,w)dTdujdwdvdx (3.8) 

[V — WjljJ 



and 



^4m)= / f*{t,x,v) 

JR»xR3xR3x52 xR+ 

xQ(f)'^(t, X + (v — w)t =f au! + T—^ — ^ — ,w)dTdu!dwdvdx. (3.9) 

[V — w)oj 

By integrating (|3.5j) with respect to variables (r, lj, v, x), we can obtain the fohowing identity 

^^^^p^ = -d^[fm+3tm)+nm) (3.10) 

where is a functional as given by 1)3.2(1 and Jf (i = 1, 2) are two functionals defined by (|3.8|1 
and ()3.9|1 . We also have a similar identity 

d^-[f]{t) 



dt 



-d-[f]{t)+3^m)+3^[m) (3.11) 



where ^ is a functional as given by (|3.2() and (i = 1, 2) are two functionals defined by (|3.8|1 
and dnSl). Estimation of (|^ and (|^ gives 

+ < 2||/o||Lir[/](i) (3.12) 

since (5(/)*(t, x, f ) < (5+(/)*(t, x, f ). Combining (|3.1U|) . (|3.11|) and (|3.12|) . we hence know that 
(j3.4j) holds. This completes the proof of this theorem. □ 

By Theorem 13.21 we then get 

Theorem 3.3. Put foo{x,v) = fo{x,v) + f^^Q{f)'^{s,x,v)ds. Assume that f = f{t,x,v) is 
a nonnegative classical solution to the Enskog equation with il.S^) and il.,'^) through an 

initial datum /o = fo{x,v) satisfying ||/o||li < 1/4 and that f € L^(R'^ x R^). Then f{t,x,v) 
converges in L^(R^ x R^) to foo{x — vt,v) as t ^ 0. 

Proof. By Theorem O we know that + (1 - 4||/o||Li)5'[/](t) < where D[f]{t) and 

5'[/](t) are the same as in Lemma 13.21 It follows that 

S)[/](t) + (1 - 4||/o||^i) fmis)ds < S)[/](0) 

Jo 
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for any t S R+, which imphes that J-^^^-^i^-^^Q^{f){t,x,v)dxdvdt < +00. Hence it can be 
easily found that 

/• + 00 p 

\f'^{t,x,v) — foQ{x,v)\dxdv< / / [Q~^{f){s,x,v) + Q~ {f){s, x,v)]dxdvds ^ 

as t ^ +00. Thus f{t,x,v) foo{x — vt,v) in L"^ as i — > +00. This completes our proof of 
Theorem □ 

In Theorem the time-asymptotic behaviour of any solution to the Enskog equation is in 
fact the time-asymptotic convergence of this solution in the norm to the free motion as t 
trends to infinity. The time-asymptotic convergence in the L°° norm is shown by Polewczak [2]. 

4. Stability 

In this section some new functionals are constructed for the stability of global classical 
solutions to the Enskog equation (|1.H) with (|1.2j) and (|1.3j) . One of them is a Lyapunov functional 
and it is equivalent to the distance functional. The time-decay property of the Lyapunov 
functional is also shown for the stability. 

Let us begin by constructing two functionals il[f,g]{t) and ^d[f,g]{'t) as follows. £[f,g]{t) is 
denoted by 

£[/, 9] (t) = { 1 + h (2) [/] (t) + 3 [g] m [/, g] (t) + /caS)^ [/, g] (t) (4.1) 

where ki and k2 are positive constants to be determined later, D is the same as given by (|3.1|l . 
£d is denoted by 

£d[f,gm= [ \f-g\*{t,x,v)dxdv (4.2) 

JR3xR3 

and S)d [/, g] {t) = 2)+ [/, g] (t) + 2)^ [/, g] (t) with 

^dif^am = [ \f-g\*{t,x,v) 

JR3xR3 



/ {f + g)"^ {t, X + t{v — w) ^ auj + t- — , w)dTdujdw 

J-R?x-R.+ xSl {v-w)u} 

:Sd[f,9m is defined by Uf,9]{t) = ^1t[f , + ^~d[f , 9m with 



dxdv. (4.3) 



^d[/'9](t)= / {v-w)u;\f-g\*{t,x,v) 



R3xR3xR3x55. 



x(/ + g)*{t, x + t{v — w)^ auj, w)dujdwdxdv . (4.4) 

This functional £ is here called Lyapunov functional. Then we have the following property of 
the equivalence between the Lyapunov and L} distance functionals £ and 2,d- 

Lemma 4.1. Let £ and £,d be defined by | |^.i| ) and ^4-^ respectively. Assume that f = f{t, x, v) 
and g = g{t, x, v) are two nonnegative classical solutions to the Enskog equation with 
and / I j.^) through initial data Jq = fQ{x,v) and go = gQ{x,v) respectively, and that f and g are 
elements in L^(R'^ x R^). Then 



uf,gm < m,9m < coUf,9m (4.5) 

2 I l|2 
Li + \\90\\li 



for any t G R+, where Cq = 1 + 2ki{\\fo\\j, + llsollri + 2A;2(||/o||li + IIsoIIlO- 
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Proof. The definition H4.1() of £ can be rewritten as 

-C[/,ff](i)= / \f -g\*{t,x,v)W*{t,x,v)dxdv 

where W{t, x, v) is defined by 

W{t,x,v) = l + kim]{t)+^[g]{t)] 

f V — w 

+k2 / {f + g){t,x — aio + T- — ,w)dTdu!dw 

yR3xR+x52 {v-w)uj 

f V — w 

+k2 I (f + g)(t,x + auj + T- r — ,w)dTdujdw, 

Jr^xK+xsI {v-w)uj 

and by Lemma l3?T] it can be found that W{t,x,v) is bounded by the norms of /o and go as 
fohows: 

1 < Wit,x,v) < 1 + 2fci(||/o||ii + IbolliO/a' + 2A:2(||/o||li + ||<7o||li)- 
Our proof of this lemma hence ends up. □ 

For any solutions / = f(t,x,v) and g = g{t,x,v) to the Enskog equation 1)1.11) with ()1.2|) 
and (|1.3|) . the time-decay properties of the two functionals D(i[f,g]{t) and can be also 

obtained as follows. 

Lemma 4.2. Let £, and £,d be defined by i4.1^ and \4-^ respectively. Assume that f = f(t, x, v) 
and g = g{t, x, v) are two nonnegative classical solutions to the Enskog equation / li. j|) with 
and M.S\) through initial data Jq = fo{x,v) and go = gQ{x,v) respectively, and that f and g are 
elements in L^(R^ x R'^). Then 

d2a[f,gm ^ 2 



dt 

and 

d^d[f,gm 



<a'Uf,9m (4.6) 



+ [1 - 2(ii/oiiii + \\go\\L^)]Uf,9m < 2{m{t)+mmuf,9]{t) (4.7) 

for any t G R+, where ^ and are given by iH.^) and i4-4\ ) respectively. 
Proof. Let us first prove ()4.6)) . By (|2.16)) . we have the following inequalities 

and 

"''*"'"''"g/*"-"''"' <m[/.#(t...„) (4.9) 
for any nonnegative solutions / and g to the Enskog equation, where 9^ is denoted by 

5^[/,5](i,a;,t>) = — / {v - w)uj[\f - g\{t,x,v'){f + g){t,x - auj,w') 

+ 1/ - 9\{t,x - auj,w'){f + g){t,x,v') 

+\f - 9\{t,x,v){f + g){t,x + auj,w) + \f - g\{t,x + auj,w){f + g){t,x,v)]dujdw. (4.10) 
Using and we have 

dUf,9m 



dt 



< I m[f,g]*{t,x,v)dxdv = a'dd[f,g]it). 



R3xR3 
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Hence 1)4. 6|) holds. Then we below prove ()4.7|) . To do this, we have to calculate 
We first calculate ^^i^i^M. To do this, we have to deduce that 



dt 



if - g)*{t, X, v){f + g)*{t, x + {v-w)t-au: + T 



V — w 
(v — w)uj 



,w) 



< Or 



(v — w)uj{f — g)'^{t, X, v){f + g)*{t^ x + {v — w)t — auj + t 



V — w 



(v — w)uj 



,w) 



+m[f, g]*{t, X, v){f + g)*{t, x+{v-w)t-auj + T 



V — w 



(v — w)uj 



,w) 



+ (/ - g)*{t, X, v)[Q{f) + Q{g)]*{t, x + [v - w)t - aoj + r 



V — w 



[v — w)uj 



,w) 



(4.11) 



and 



(g - f)*{t, X, v){f + g)*{t, x + {v-w)t-aoj + T 



V — w 
(v — w)uj 



,w) 



< Or 



{v - w)uj{g - f)*{t, X, v){f + g)*{t, x + {v - w)t - au + ^ 



V — w 



(v — w)uj 



,w) 



+m[f, g]*{t, X, v){f + g)*{t, x+{v-w)t-auj + T 



V — w 



(v — w)uj 



,w) 



{g - f )*{t, X, v)[Q{f ) + Q{g)]*{t, x + {v - w)t - auj + r 



#1 



V — w 



-,w) 



(4.12) 



{v — w)uj ' 

for any fixed variables {x,v,w,uj). In fact, (|1TTT|) and (|1TT^ result from ^J^, (ESJ, (|S?7)) . 
(|4.9() and the two identities as follows: 



dt 



g'^it, X + [v — w)t — aui + T 



V — w 
{v — w)uj 



w 



{v — w)dxg'^{t, X + {v — w)t — au) + T 



V — w 



{v — w)uj 



w) 



+Q{g)*{t, X + {v — w)t — auj + T 



V — w 
(v — w)uj 



w) 



and 



(v — w)dxg'^{t, X + {v — w)t — auj + T 



V — w 
{v — w)uj 



w 



= {v — w)ujdr 

Denote Zf[f,gm (i = l,2) by 



g'^it, X -\- {v — w)t — au + T 



V — w 



R3xR3xR3x5ixR4 



{v — w)oj ' 
^[f,g\*{t,x,v) 



Wj 



X (/ + g)*{t, x + {v — w)t =F aw + r 



V — w 



, w)dTdijjdwdvdx 



and 



R3xR3xR3xS'ixR4 



{v — w)u 

\f -g\*{t,x,v) 



(4.13) 



X (Q(/) + Q{9))* (t, x + {v-w)tTauj + T 



V — w 
{v — w)oj 



w)dTdujdwdvdx. (4-14) 
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Then, using H4.11() and H4.12() . we can obtain the following inequality 

< -St if, 9m+ if, 9m+ [f,9m. (4.15) 



dt 

Similarly calculating ff^i-^i^ll^l ^ -^^^g also have 



'^^''^l^f^^^ < -r, [/, 5] ) + [/, 9m + ^2 [f,9m. (4. le) 

Combining (|4.15j) and (|4.16j) . we get 

"^^'^l^;^^^^^ < -uf,9]{t) +3r[/,5]w +^2[f,9m. (4.17) 

By (|4.13j) and (|4.14jl . the estimates of 2t[fi9]{'t) + ■S^[f^9]{'t) {i = 1)2) can be obtained as 
follows: 

^t[f.9m+z^[f,9m < 2(ii/o||li + \\9o\\L^)si[f.gm (4.18) 

and 

2t[f,9]{t) +Z2[f,9m < 2mm)+d[9Kmd[f,9m. (4.19) 

By applying (|4.18() and ()4.19() into the estimation of the right side of H4.17() . (|4.7() follows. □ 

It can be further shown that the functional £, has the following time-decay property. 

Theorem 4.3. Let 2 and be defined by i4-l\ ) o,f^d \4-4\) respectively. Assume that f = 
f{t,x,v) and g = g{t,x,v) are two nonnegative classical solutions to the Enskog equation 
with m.l^) and M.clp through initial data fo = fo{x,v) and go = gQ{x,v) satisfying ||/o||li < 
(2 — \/2)/4 and II (7o II Li < (2 — \/2)/4 respectively, and that f and g are elements in (R? x H?) . 
Then there exists a positive constant a G (0, 2-v/2) such that 



d£[f,g]{t) . , Q 



dt 



2 



4(2 + a) 



k + 1 



Uf,9m (4.20) 



with ki = (2-\-a)k2 and k2 = ka^ for any t € R+, where k > -^z^- Furthermore, the Lyapunov 
type estimate holds as follows: 



£.[f,g]{t) + C f ^d[f,g]{s)ds<il[f,gm (4.21) 
Jo 

for any t G R+, where C is a positive constant independent oft. 

Proof. By the definition 1)4. 1|) of the functional £, using Theorem 13. 21 and Lemma 14.21 we have 

^^^M(!) = [1 + ,,(^[/](t) + mm '-^^^^ 

< [a\i{m\{t) + mm) + 2A:2(||/o||li + Ibollii) - ^2 + 0^] ^d[f,9m 

+ [Aki\\fo\\L^-ki+2k2]mmd[f,9]{t) + [ikM\Li-ki+2k2]mmd[f,9]{t)- 

Since S)[/](t) < ^ll/ollii and D[g](t) < ^Ibollii, it follows that 
d2.[f,g]{t) 



< 



[2A;i(||/o||ii + Ibollii) + 2A:2(||/o||li + IbollLO - ^2 + a'] Uf,9m 



dt 

+ [4A:i||/o||li -h + 2k2]d[f]md[f,9]{t) + [4A;i|bo||Li -h + 2^2] %](t)i2d[/, sK*)- 
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Notice that ||/o||li < (2 — \/2)/4 and Hs'oIIli < (2 — V^)/4. There exists a positive constant 
a G (0,2\/2) such that ||/o||li < i(JXa) IbollLi ^ 4{2+«y' ■'^^^ choose /ci = (2 + a)A;2 and 

k2 = ka? where k > ^3^. By estimating the right side of the above inequahty, it can be found 
that 

" ^ ,2 



d£[f,g]{t) . , a 



dt 



< a" 



2 



+ 1 



_4(2 + a) 

(|4.20|) and (|4.2H) then follows. The proof of Theorem 14.81 is hence finished. □ 

By Theorem 14.81 we further have 

Theorem 4.4. Assume that f = f(t,x,v) and g = g{t,x,v) are two nonnegative classical 
solutions to the Enskog equation U.l\) with MJJ\) and M.'J\) through initial data /o = fo{x,v) and 
90 = go{x,v) satisfying ||/o||li < (2 — V^)/4 and H^oHli < (2 — respectively, and that f 

and g are elements in L^(R'^ x R^). Then 

[ \f-g\it,x,v)dxdv<C [ \fo-go\ix,v)dxdv (4.22) 

JR3xR3 JR3xR» 
for any t G R_|-, where C is a positive constant independent oft. 

Proof. By Theorem 14.81 there exists a positive constant a G (0,2^/2) such that £[f,g]{t) < 
£[/,ff](0) with ki = {2 + a)k2 and k2 = ka"^, where k > By using Lemma l4. 11 it follows 

that 

Uf,9m < ^[f,gm < £[/,<?](0) < Co£d[/,5](0). 
This hence completes the proof of Theorem 14.41 □ 

Theorem 14.41 can be here regarded as an application of the Lyapunov functional £ into re- 
covering a result of Cercignani about the stability of classical solutions to the Enskog 
equation under the assumption of the suitably small initial data. 
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Appendix A. A Counterexample of Inequalities Given by Ha 
Below is shown a counterexample of two inequalities as follows: 

/ /"^{x + t{v — v^) ± aa + Tn{v,v^,),v^,t)dTda < / f'^{x,v^,t)dx (A.l) 

Jslx-R+ Jrs 

where S\ = {a : (u — f*)cj > 0, o" G 5^}, n{v, v^) = {v — v^:)/\v — v^\ and f{x, v, t) is a nonnegative 
function of the space x G R^, the velocity v G R^ and the time t G R+. These inequalities are 
given by Ha (see (3.1) and (3.2) at Page 1004 in |^) and applied into the proof of his theorems 
which show that two so-called Lyapunov functionals have their Lyapunov type estimates that 
yield the time-asymptotic behaviour and the stability of global classical solutions to the 
Enskog equation without external forces. 

To give a counterexample, we first have to prove the following lemma. 

Lemma A.l. Let F{x) be an integrable function on R^ and n a unit vector in R^. Assume that 
S\ = {oj : nuj > 0,a; G S^} where S"^ is a unit sphere surface in R'^. Then 

I F(x)dx = F{auj + Tn)nujdTduJ (A.2) 

where a / 0. 

Proof. Without loss of generality, we assume that a > 0. Let u) = (sin cos sin sin 99, cos 0), 
where 9 is an angle between v and u:,0<6<^,0<ip< 2tt. Denote n and x by n = (ni, n2, na) 
and X = {xi,X2,xs) respectively. Take a transformation x = au + rn, i.e., 

xi = TUi + asinO costp, X2 = Tn2 + as'mO s'mip, x^ = rn^ + a cos 9. (A. 3) 

By direct calculation, it then follows that the Jacobian determinant of this transformation is 
s'm9. HA.2|) hence holds. Our proof of this lemma is thus finished. □ 



'{xi,X2,X3) 



a nuj 



d(r,e,ip) 

Lemma IA. II not only implies that F{x) must be integrable in R'^ if F{auj + rn) is integrable 
in R X S'^, but also indicates that F{auj + rn) might not be integrable in R x S\ even if F{x) 
is integrable in R^ since it can be known from HA.3|) that nuj = {xn — T)/a and (r — xn)^ = 



(xn)^ + a? — x^ . For example, let us choose n = (1,0,0) and F(x) = F(x) j \J\x\ + Xg — a^l, 
where F(x) G C^(r2), F(x) = 1 as x G i?2a5 < F(x) < 1 as x G ^2, B^a C 0, ^2 is a compact 
subset in R^ and i?2a represents a ball with its center in the origin and its radius 2a. It is easy 
to see that Fiauj -|- rn) is not integrable in R x S\ but F(x) is integrable in R^. This is a good 
counterexample of inequalities (|A.1|) and so the two inequalities do not hold mathematically. 
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